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Wo characterise the primordial perturbations produced due to both inflaton and curvaton fluc- 
tuations in models where the curvaton has a quadratic, cosine or hyperbolic potential, and the 
inflaton potential is characterised by the usual slow-roll parameters. Isocurvature curvaton fleld 
perturbations can produce significant non-Gaussianity in the primordial density field, in contrast 
with adiabatic infiaton field perturbations which produce negligible non-Gaussianity for canonical 
scalar fields. A non-self-interacting curvaton with quadratic potential produces a local-type non- 
Gaussianity that is well described by the non-linearity parameter /nl, which may be scale-dependent 
when the inflaton perturbations dominate the power spectrum. We show how observational bounds 
on non-linearity parameters and the tensor-scalar ratio can be used to constrain curvaton and infla- 
ton parameters. We find a consistency relation between the bispectrum and trispectrum parameters 
in a mixed inflaton-curvaton model for a quadratic curvaton potential. Self-interaction terms in the 
curvaton potential can lead to both a large trispectrum parameter, ^nl, and scale-dependence of 
the non-linearity parameters. 

I. INTRODUCTION 

Inflation is our most successful theory for explaining the initial conditions required for the hot Big Bang cosmology. 
In particular, primordial density perturbations can be produced from initial quantum fluctuations that are stretched 

by the accelerated expansion up to supcr-Hubblc scales to become the large-scale structure of the Universe today. 
This mechanism can give rise to an almost scale-invariant power spectrum, as observed in the cosmic microwave 
background [1]. 

The standard model of inflation is described by a light scalar field 0, the inflaton, whose slow-roll controls the 
potential energy that drives the accelerated expansion. When inflation ends the inflaton potential energy rapidly 
decays into radiation. Any light scalar field during slow-roll inflation will acquire a nearly scale-invariant spectrum 
of perturbations at Hubble exit, k = aH, and any such field can, in principle, source structure in the Universe. The 
curvaton, is a light, weakly- coupled field during inflation whose energy density is negligible during inflation, but 
if the field remains weakly coupled at the end of inflation, its energy density can grow relative to radiation after 
inflation, and perturbations in the curvaton fleld can lead to primordial density perturbations when the curvaton 
decays into radiation [2-8]. A distinctive feature of the curvaton model is that it can leave behind signiflcant local- 
type non-Gaussianity in the primordial density distribution, |/nl| ^ 1- 

Previous studies [9, 10] have used primordial non-Gaussianity, characterised by the non-linearity parameter /nl, 
and the tensor-to-scalar ratio, rr, to constrain curvaton model parameters. In the simplest model of a curvaton with 
a quadratic potential, the amplitude of primordial density perturbations, together with a measurement of /nl and 
would fix the energy scale of infiation, the initial curvaton vacuum expectation value (VEV) and the dimensionless 
curvaton decay rate. A self-interacting curvaton would introduce additional model parameters that could be fixed 
by scale-dependent /nl and higher-order correlators. These studies focussed solely on density perturbations coming 
from the curvaton field. But if inflation is driven by a slow-roll inflaton, then there are inevitably fluctuations in 
the inflaton field too which lead to some level of density perturbations when inflation ends and the inflaton energy 
decays into radiation. Inflaton perturbations are adiabatic and thus lead to constant curvature perturbations, on 
super-Hubble scales, whose local-type non-Gaussianity is always small, |/nl| < 

In this paper we also consider the possibility that both fields contribute for the primordial density perturbation. 
In particular we note the possibility that while the infiaton contribution to the primordial power spectrum may 
dominate over that from the curvaton, the curvaton can nonetheless source significant non-Gaussianity. Bartolo 
and Liddle [11] were the first to consider in what regime the linear curvaton or inflaton fleld perturbations would 
dominate the primordial power spectrum in a simple model of two massive scalar fields. Ichikawa et al [12] also 
considered non-Gaussianity and the tensor-scalar ratio in several classes of inflation models, including chaotic, hybrid, 
and new inflation. In this paper wo will evaluate the relative contribution of curvaton field fluctuations to the 
primordial power spectrum for a general inflaton potential. The contribtition of adiabatic inflaton perturbations is 
given relative to the tensor power spectrum by the slow-roll parameter, Here we define an analogous parameter, 
Cc, describing the contribution of curvaton perturbations relative to the tensor power spectrum. The second inflaton 
slow-roll parameter, r], then affects only the scale-dependence of the power spectrum. Very recently Kobayashi and 
Takahashi have considered the scale-dependence of both the power spectrum and the non-Gaussianity including mixed 
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inflaton-curvaton models [13]. Note that Langlois and Vernizzi [14] studied linear adiabatic and isocurvature density 
perturbations allowing for a general inflaton potential, which was then extended [15] to second- and higher-order 
to study non-Gaussianity in both adiabatic and residual isocurvature density perturbations after inflation and their 
correlation. Recently Kinney et al [16] considered the complementary constraints coming from the shape of a mixed 
(bimodal) power spectrum for adiabatic density perturbations. In the following we assume that all species are in 
thermal equilibrium after the curvaton decays and thus there are no residual isocurvature modes [17, 18]. 

In Section II we will briefly review the origin of density perturbations coming from both inflaton and curvaton 
field perturbations, and the spectrum of tensor metric perturbations (gravitational waves). We will review both 
the linear transfer of field perturbations into radiation, and the non-linear transfer for curvaton perturbations at 
second- and third-order in the field perturbations which can give rise to a non-vanishing primordial bispectrum and 
trispectrum. In Section III we present our results based on numerical solutions of the curvaton field evolution after 
inflation [10] and previous numerical studies of curvaton decay based on a fluid description of the curvaton at late 
times [19 22]. We focus primarily on the simplest curvaton model with a quadratic potential with a fixed mass, m^. 
In this case signiflcant non-Gaussianity arises only when the curvaton is sub-dominant when it decays, ^ 1 and 
we derive a consistency relation between the bispectrum and trispectrum parameters which holds even in the mixed 
inflaton+curvaton case. The third-order non-linearity parameter, ^tnl, remains small, ^nl /nl^, even in the mixed 
scenarios. Scale-dependence of /nl may distinguish the mixed inflaton-|-curvaton model from the curvaton limit for a 
quadratic curvaton potential. We also examine self-interacting curvaton models, including a cosine potential, which 
introduce an additional mass scale, /, where self-interaction terms become important for x* ~ /• Self-interactions 
can produce significant non-Gaussianity even if the curvaton dominates when it decays, ~ 1. Self-interacting 
curvatons produce large third-order non-linearity parameter, ^nl, as well as scale-dependent /nl- We discuss our 
results and conclude in Section IV. 



II. DENSITY PERTURBATIONS FROM INFLATION 



A. Background evolution 

We take both the infiaton and the curvaton to be in a slow-roll regime during inflation but assume that the 
Friedmann equation is dominated by the infiaton potential energy 
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where mpi is the reduced Planck mass, rupi = SttGn 
We define the slow-roll parameters 



1 2 fVA 



and 

— 2 ^AB ... 

Vab = mpi—, (4) 

where Va = dV/dA. Note that in the slow roll approximation e ~ <C 1 and in the curvaton scenario we assume 

that <C 60 so that e ~ e^. We also assume the fields are decoupled so that r]^^ = 0. 



B. Perturbations during inflation 

During inflation, any light scalar fields (with effective mass less than the Hubble scale, |?7| < 1) acquire a spectrum 
of perturbations due to vacuum fiuctuations on sub-Hubble scales being stretched up to super-Hubble scales by the 
accelerated expansion. In particular the curvaton and inflaton fleld perturbations on spatially-flat hypersurfaces at 
Hubble exit have a power spectrum 

~ Pax* - (^) ' 
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where we neglect slow-roll corrections, including the cross-correlation between inflaton and curvaton perturbations 
[23]. 

Since the inflaton determines the energy density during inflation, inflaton field perturbations on spatiaUy-flat 
hypersurfaces, 5(j), correspond to adiabatic curvature perturbations on uniform-density hypersurfaces at Hubble exit, 
C* = —{H54>/(j))^ and hence we have 

On the other hand curvaton fluctuations are isocurvature fleld perturbations during inflation, <C e^, and remain 
eff'ectively frozen, x — 0, and hence are gauge-invariant during inflation. In particular we can identify curvaton fleld 
perturbations on spatially flat hypersurfaces with relative entropy perturbations [24] 

Sx* 0^ {^X - |<50^ - 5x* ■ (7) 

The power spectrum of free gravitational waves (tensor metric perturbations), like light scalar flelds, only depends 
on the inflation scale at Hubble exit 

^- = ^ff^)'. (8) 

The tensor-to-scalar ratio is defined by 

- ^ (9) 

Both the adiabatic curvature perturbation and the tensor perturbations remain constant on super-Hubble scales, so 
we have a tensor-scalar ratio during inflation [25] 



Vt 



= 16e* . (10) 



The tensor spectral index is due solely to the variation of the Hubble scale during inflation 

However the inflaton field and curvaton field evolve on super-Hubble scales due to their effective mass, and gravitational 
coupling for the inflaton field, so their spectral tilts are given to leading order by [25] 

n<p-l^ ^=-6e*+2^,, (12) 



C. End of inflation and after 



At the end of inflation the inflaton decays completely into radiation transferring its curvature perturbation to the 
radiation, = C^. We assume that reheating or preheating does not alter the power spectrum of the adiabatic density 
perturbation on large scales, nor does it alter the fluctuations of the curvaton fleld on large (super-Hubble) scales. 

The curvaton stays in an over-damped regime until the Hubble rate drops to ~ m^. At this point the curvaton 
starts oscillating and behaves like a prcssureless matter fluid. (We will not consider the possibility of the curvaton 
driving a second period of inflation [26], i.e., we assume x* < wp;.) Once the curvaton starts evolving like a 
pressureless fluid we can write its local energy density on uniform-total-density hypersurfaces, p^, in terms of its 
homogeneous value, p-^, and the inhomogeneous entropy perturbation [15, 27] 



Px = Pxe^^^""^*' = Pxe^" 



(14) 
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where = 3 (Cx ~ C(p) is the non-adiabatic part of the curvaton perturbation. One should note that in the standard 
curvaton scenario one takes ^ C^, hence S-^ ~ 3^;^. In the mixed inflaton-curvaton case this may no longer hold, 
therefore the quantity to use is S-^ ■ When the expansion rate drops to iJ ~ the curvaton starts oscillating in the 
bottom of its potential, behaving like a pressureless, non-interacting fluid. At later times, but before the curvaton 
decays, the potential of the curvaton field can be well approximated by a quadratic potential and its time-averaged 
energy density can be described by 

,2 I, ,2 



Px = ^K\x\- (15) 



One can use Eq. (14) to determine the relation between the entropy perturbations of the curvaton and its field 
fluctuations during inflation. In the beginning of oscillation we have 

Pxe^" = Irnlxlsc (16) 

where — "Ti^Xosc/^- Note that the subscript osc stands for beginning of oscillations. Let's define Xosc = 9ix*) 
where g accounts for non-linear evolution of the field between inflation and oscillations [22] . If the curvaton potential 
is quadratic and we can neglect the self-gravity of the curvaton, we expect linear evolution. On the other hand if it 
is not quadratic throughout all evolution we need to correct the fleld perturbations. It is convenient to expand Xosc 
in terms of field perturbations during infiation, 5x*, i-e., 

Xosc -9 + 9'5X* + \9"5xl + ■ • ■ , (17) 

where primes denote derivatives with respect to x*- Expanding both sides of Eq. (16) up to second order we find that 
the curvaton entropy perturbation is 



<5x = 2^^X* 



9" / 9 



5xl + 0{5xl). (18) 



One should note that for a non-interacting, isocurvature field, 5x* is a Gaussian random field. Therefore we can 
separate the curvaton entropy perturbation in a Gaussian linear part and second order term as 

Sx=SG + -A'4T-l\sh (19) 



4 Vfl" 
where 

Sg = 2^-5x* . (20) 

Hence, using Eqs. (5) and (20), the power spectrum of entropy perturbations in the curvaton is, at leading order, 
given by 

9'V (H,-' 



^s.=^So=H'-} ■ (21) 



D. Transfer of linear perturbations 

The curvaton decays into radiation when ~ F^. We will consider that the curvaton (and the inflaton) decay 
prior to CDM freeze-out. Therefore we won't consider any residual isocurvature perturbations after curvaton decay 
[17]. 

The primordial density perturbation produced by curvaton decay can be estimated using the sudden decay ap- 
proximation [5]. This assiimcs that the curvaton happens instantaneously on the total- uniform-density hypersurface 
H = Tx- Before the curvaton decays, C7 = C^- Therefore we know that on the sudden decay hypersurface we have 



(22) 
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At sudden decay we have that the final radiation energy density p = + p-)^. Hence after decay we have 

(1 - Ox) e^(f*-^) + O^e^ff^-f) = 1 . (23) 

After the decay we have a constant curvature perturbation on super-Hubble scales. Expanding Eq. (23) to first 
order, we have 



C = i«xCx + (l-«^- (24) 

where [5, 8] 



C^ + ^^x- (25) 



(26) 



dec 



Since the adiabatic infiaton field perturbations and the isocurvature curvaton field fluctuations (5) are uncorrelated, 
the power spectrum of the total primordial curvature perturbations, after curvaton decay, is given by 

Vc=Vc, + ^Vs,, (27) 
Following Eq. (6), and using Eq. (21), we can write this as 

1 /I 1 \ [H^\ 

^^ = 2. 

where we define a quantity 



8 \g'mpi J R. 



' 1 



2 • (29) 



■-x 

The curvaton contribution to the primordial power spectrum corresponds to (2TOp;ec)~'^(-ff*/27r)^, i.e., plays the 
same role for the curvaton contribution to the final power spectrum as e« does for V^^ in Eq. (6). Thus Cc marks the 
critical value of e* between infiaton-domination of the primordial power spectrum and curvaton-domination of the 
power spectrum. It follows that we can write 

^ ■ (30) 

fc* ~r fcc 

can be seen as the function that weighs the curvaton contribution to the final power spectrum. For e* Cc 
the curvaton is the dominant contributor to scalar perturbations and w-^^ ~ 1. In the opposite regime, e* <C Cc, the 
infiaton dominates the primordial power spectrum and <C 1. 

The spectral index of the primordial power spectrum is then given by 

nc-l= = w^{n^ -!) + (!- w^){n^ - 1) (31) 

Substituting the tilts (12) and (13) for each field in (31) we have [14, 25] 

- 1 = -2e. + 2r]^^w^ + (1 - w^){-4e, + 2r,^^) . (32) 

The running of the power spectrum, assuming slow roll infiation and neglecting curvaton-infiaton interactions, is 
given by 

= = ^«x"x + (1 - w'x)"0 + w'x (1 - w'x) (% - n<t'f (33) 

with [28, 29] 

= 166*7700 - 24e2 - 2e' (34) 
= 4e* (-2e» + -q^^ + rj^^) , (35) 

where & = {d'^V/d(p'^)/9H*. 
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E. Non-linearities 



Using the SN formalism we identify the non-linear curvature perturbation, ^, with the perturbed expansion up to 
a final uniform-density hypersurface, A'', as a function of the local field values on super- Hubble scales during inflation 

[30] 



(36) 



A,B 



A,B,C 



where = dN/dcf)^ and N^b = d? N / dcf) j^dcf) b , etc. We define the first non-Gaussianity parameter in terms of the 
amplitude of quadratic contribution to C, relative to the linear terms: 



/i 



NL 



5 Y.abNaNbNab 

6 NaNbSab? ' 

If we consider terms in Eq. (23) up to second order we find [15] 



c 



Plugging (19) into (39) we find 



/■ /- ^X o ^X 



18 



1 



99 

„i2 



It is straightforward to see from Eqs. (6), (21) and (27) that 

25' 



R 



1 



(37) 

(38) 
(39) 

(40) 

(41) 
(42) 



We only need to consider the linear terms from the inflaton since 7V^^ <C (-/V^)^ and N^j,^ = 0. Then, the first 
non-Gaussian parameter (37) for curvaton-|-inflaton simplifies to [15] 



f _ 5 ^XX „„2 



It follows directly from Eq. (40) that /nl is given by 



/NL 



99 

g,2 



2, & ^ 



w: 



X ■ 



Note that taking the derivative of Eq. (41) with respect to x* we get 



-'^XX ~ g^X 



9 5^ 



R 



Then, using Eq. (43) we find the general expression for /nl in the mixed curvaton-infiaton scenario 



/nl = 



4i?^ 



9 9 

„/2 



+ 



5 R^'g/g' - 2R^ 
4 i? 2 



Comparing this with Eq. (44) obtained in the sudden-decay case we have [22] 



JrC-y — - — ZJrC-v — —JrC-v — —tt-y 

g' 



(43) 



(44) 



(45) 



(46) 



(47) 
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The third order non-Unear parameters are gfNL and tnl- They are defined by [30, 31] 

^ 2b Y.ABcNANBNcNABC 

E^BC NaNbNacNbc 
[HabNaNbSabT 

For the inflaton+curvaton case Eqs. (48) and (49) reduce to 



5nl 



54 iV3 



7V2 

Tnl = -^w. 



7V4 '"x • 

X 



Taking the third derivative of A'' with respect to x* we find 



9 9^ 9 



2 X,dec 

Substituting Eqs. (41) and (52) into Eq. (50) we get 

1 



5nl 



25 
24 



V3 



5 

5'2 



+ 2 +2 



n'3 



+ 



Rx'9'\ 



Using Eq. (47) for the sudden-decay approximation we can ehminate the derivatives of to obtain [22] 



25 
54 



4i?. 



_£ L o p 9 

'3 ^^g>2 
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1 + 



+ 



5'- 



lOi?;^ + 3i?x' 



(48) 
(49) 

(50) 
(51) 

(52) 

(53) 

(54) 



Equations (46) and (53) do not rely on the sudden decay approximation. Nonetheless Ref. [22] showed that the 
sudden-decay formulas (44) and (54) do give a good fit to /nl(^x) ^^'^ 5nl(-Rx) from the full numerical solution 
with continuous decay. For example, /nl(.Rx) accurate to 1% for /nl > 60. Therefore in the following we will use 
Eq. (44) and (54) to give the non-linearity parameters as a function of R^. 

Note that from Eqs. (43) and (51) we have 



g/NL 



(55) 



The inequality tnl > (6/nl/5)^ is an important test of non-Gaussianity the mixed curvaton+inflaton scenario and 
multi- field scenarios in general [35] with equality only in the curvaton limit, w-^^\. 



III. NUMERICAL RESULTS 



Wo now wish to compute observable quantities such as /nl, <7nl and tt for different model parameters. We assume 
that the curvaton is efi'ectively frozen during infiation, which should be a good approximation while the effective 
curvaton mass is much less than the inflationary Hubble scale, i.e, 77^^ <C 1. We will then numerically solve for the 
local evolution of the curvaton field during the radiation- and, possibly, curvaton-dominatcd epochs after inflation, 
until the curvaton starts oscillating in the minimum of its potential and behaves like a pressureless fluid, but before 
it decays. We allow approximately 10^ oscillations, i.e., we assume sufficiently slow decay, F/m < 10"^^ consistent 
with the hypothesis that the curvaton is weakly coupled to other fields. 

We numerically solve the Klein-Gordon equation for x prior to decay 

X + 31fx + ^x = 0, (56) 

where the Friedmann equation takes the form 



(57) 
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and the curvaton density is given by = + V{x)- This allow us to take different potentials for the curvaton 
field. 

The initial conditions to solve Eq. (56) are Xi — X* ^^id Xi — ^V^/^Hi, since the curvaton is slow-rolling down 
its potential. We take the universe to be radiation dominated initially, p-y^i ^ p^^i and hence Hf 3> T^(x*)/3mp;. In 
practice we set p-y^i = lO^'^m^GeV^ in our numerical solutions, where m-,^ is the mass of the curvaton at late times. 
This ensures that Hi > lOOm-^ and p^^i m^Xi for any Xi < mpi- 

Once the curvaton starts oscillating we can compute the quantity [10] 

which becomes constant during oscillations as — )• const/a^ and we have = Xosc <^ Pfw- In this way we connect 

the scalar field description of the curvaton with a fluid description which has previously been used to numerically 
study the decay of the curvaton [19-22]. Following [10] we can compute the efficiency parameter i?^ using the fitting 
formula [20] 

^ . 0.924 Prru: 

R^ippw) Y^Jy^PpwJ . (59) 

We find that, for a given curvaton potential, ppw is dependent only of the initial curvaton field value, Xi- Therefore 
is dependent only upon the initial curvaton field value and the dimensionless decay rate, T^/m^. 

Going beyond our previous work [10] we will include inflaton perturbations in addition to curvaton field pertur- 
bations in our computation of the primordial density perturbation. However the inflaton perturbations represent 
adiabatic perturbations on super- Hubble scales, i.e., local perturbations along the same background trajectory [24], 
and they can be treated independently of the curvaton field perturbations. Looking at the total power spectrum, 
Eq. (27) we see that we have gained an extra free degree of freedom, e*, with respect to the purely curvaton limit 
(e, Cc). Therefore our free parameters will be T-^/m-^, x* and e*) for the quadratic curvaton [10]. Going beyond 
the quadratic curvaton potential we will consider models where self-interactions lead to the potential becoming flatter 
or steeper beyond a characteristic mass scale, /, introducing one new parameter in addition to the curvaton mass 
about the minimum of its potential. 

WMAP 7 [1] gives V(; ~ 2.43 x 10~^ for the amplitude of the power spectrum of curvature perturbations. We will 
use this observational constraint to fix the infiationary scale. Since the power spectrum (27) is proportional to 
we can find a value of if* that gives the correct power for any values of the other curvaton model parameters. Using 
Eqs. (28) and (30) we arrive to the constraint equation 

= 2V27Tmpi (e-i + e^'y'^^ V^'/^ , (60) 
where the value of Cc is determined numerically via the formula 

e.= ^-( , V- (61) 

2 \ppwmpiR^J 

The tensor-to-scalar ratio is then 

Vt = 16w;^ec = 16e*(l — w-^) . (62) 




A. Quadratic potential 

The simplest potential consistent with a curvaton scenario is a quadratic potential 

V{x) = lmlx\ (63) 

and we will see that it is also a good decription of the behavior in more general cases when the curvaton is sufficiently 
close to the minimum of its potential. 

We start by studying the relative contribution of the curvaton to the scalar power spectrum characterised by 
parameter defined in Eq. (29). In Figure 1 we plot Cc as a function of the initial curvaton VEV, %*, and the 
dimensionless decay time niy^/r-^. 
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FIG. 1: The plot show contours of ec defined in Eq. 29 as a function of the curvaton parameters, x*i the curvaton VEV, and 
m^jV^ the dimensionless decay time, for the quadratic curvaton potential, Eq. (63). Curvaton perturbations dominate the 
primordial scalar power spectrum for e, ^ ec- 



While the curvaton remains subdominant in the radiation era, we have an analytic solution for the curvaton field 
[12, 14] and we find [10] 



Vfw 



1.046x* 



(64) 



We can clearly identify the two analytic regimes in Figure 1 



1.125 f ^\ 

r 



10.. 



X* 



for i?^ ~ 1 , 
for i?^ < 1 , 



(65) 



corresponding to the straight lines in Figure 1. 

If we now include the contribution from inflaton perturbations to the total scalar perturbation, we can identify 3 
regimes of interest which depend on the value of e.^, for a given e^- 

1. The curvaton limit corresponds to Cc ^ e*. In this case most of the structure in the universe comes from the 
curvaton, i.e., ~ 1 from Eq. (30). This case has been studied in our previous work [10] and in most of the 
curvaton literature. We can identify this limit in Figure 1 for a fixed value of e* as the region inside the contours 
towards the right of the plot, i.e., for long decay times (T^ ^ "^x)- 

From Eq. (62) we have in the curvaton limit 



tt ^ 16ec 



(66) 



Therefore upper bounds on the tensor-to-scalar ratio place constraints on Cc but do not directly constrain e* 
since e* ^ Cc- 

In the curvaton limit, ~ 1, and assuming an effectively massless curvaton, t]^^ <C 1, then Eq. (32) gives a 
red spectral tilt, n,; — 1 ~ —2e^.. In this limit the tilt gives a measurement of the first slow roll parameter, e*. 
Consider a fiducial value ~ 0.96 consistent with WMAP7 [1]. For this value of ~ 0.02 we can identify the 
curvaton limit with the region to the right of the contour — 0.02 in Figure 1. 

2. The second limit of interest is Ec S> e*. This is the case for which curvaton perturbations are sub-leading in the 
scalar power spectrum, i.e., <C 1 in Eq. (30). These regions correspond to a parameter range where the decay 
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happens too fast (bottom left of the plot), or the curvaton VEV is too big (top) suppressing the curvaton power 
spectrum. The region Ec ^ 1, in Fig. 1, will always be in this inflaton dominated limit in slow-roll inflation 
since e* ^ 1. 

The tensor-scalar ratio rx directly constrains the slow- roll parameter in the this limit. From Eq. (62) we 
have 



tt ^ 16e* . (67) 

The spectral tilt of the primordial scalar power spectrum (32) is determined by the usual inflaton slow-roll 
parameters, rt^ — 1 ~ — 6e* -I- 2rj^^ for w^<^l. 

In this limit the presence of the curvaton may still be important to as a source of primordial non-Gaussianity 
or residual isocurvature perturbations after the curvaton decays [15]. 

3. The third region of parameter space corresponds to e.^ ~ Cc which corresponds to a mixed scenario. In this case 
the tensor-scalar ratio (62) no longer directly constrains e* or Cc but the combination 

rx = . (68) 

ec + e* 

For example, an observed tensor-scalar ratio, r^, places a lower bound on the slow-roll parameter, > rT/16. 

In Figure 2 we show contour plots for the non-Gaussianity parameter, /nl, and the tensor-scalar ratio, ry, for the 
case e* = 0.02. The thin black dotted line is the contour line ec — 0.02 which marks the borderline between the region 
(1) described above with curvaton-dominated primordial power spectrum and region (2), inflaton-dominated. The 
curvaton limit, region (1), lies to the right of the €c — 0.02 contour. 

We also plot the current observational upper bound on the tensor-scalar ratio, rx ^ 0.24 [1]. For a given value of 
e*, the contours of the tensor-to-scalar ratio follow the contours of Ec plotted in Figure 1, as expected from Eq. (68). 
However, rather then growing without bound as becomes large, as happens if we consider only the curvaton 
perturbations [10], in the presence of a finite e* the tensor-scalar ratio saturates with rx 16e* in region (2) where 
ec ^ e*. For e* = 0.02, for example, the tensor-scalar ratio is bounded so that rx < 0.32. 

Similarly the inflaton's (Gaussian) contribution to the primordial scalar power spectrum suppresses the non-linearity 
parameter /nl for ec > e* in region (2). We see that the largest values for /nl occur in region (3), near the 
boundary between the curvaton- and inflaton-dominated power spectra, where ec — e* . In the absence of any inflaton 
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iog,„ m/r 

^10 

FIG. 3: The plot show contour hnes for the non-Unear parameter /nl (blue lines) and the tensor-to-scalar ratio rr (thick red 
lines) as a function of the curvaton parameters, x* smd m-^/T^, for the quadratic curvaton potential, Eq. (63), and a fixed 
value of the inflation slow-roll parameter, e, = 0.1. The black broken line corresponds to tc = 0.1. 



perturbations {w^ = 1), the non-Gaussianity continues to grow without bound as x*/'^Pi — >■ for a fixed value 
of Triy^/T^ [10]. But also becomes large as x*I™'P\. ~^ ^i^d therefore the infiaton perturbations dominate the 
primordial power spectrum. From Eq. (46) we see that /nl is suppressed by an additional factor uj^ ~ £*/ec ^^"^ 
have 

/nl ^ I ^ 0.0386^ f ^ V^V-^ V , (69) 



which is suppressed as x*l'^Pi ~^ for ^ given m^/F^. 

If we demand a lower bound on the non-Gaussian parameter, /nl > 10 for example, this places an upper bound on 
the curvaton VEV, x* < 1-2 x 10^^ GeV for e» = 0.02, but also a lower bound on the decay rate F^ < 3 x 10~®m^. 

Figure 3 is similar to Figure 2 but corresponds to a larger slow-roll parameter = 0.1. The thin black dotted line 
separating the infiaton- and curvaton-dominated regions here corresponds to ec = 0.1. We see that for larger values 
of e* the parameter regime (1) corresponding to the curvaton limit extends to smaller values ofm^/T^ and larger x*, 
permitting larger values of /nl- 

On the other hand observational bounds on the tensor-scalar ratio now place more severe constraints on the allowed 
parameter values. If we put a lower bound on the non-Gaussian parameter, /nl > 10 for example, this places an 
upper bound on the curvaton VEV, x* < 2.5 x 10^^ GeV, and a lower bound on the decay rate, F^ < 2 x 10~''mj(-, 
for =0.1 and < 0.24. 

The entire infiaton dominated region (2) is excluded by observational bounds on the tensor-scalar ratio for such a 
large value of e*. On the other hand e* = 0.1 is allowed in much of the curvaton dominated region (1). 

Note however that such a large value of e* requires a similar positive value of rj^^ , tuned such that the spectral tilt 
remains small in the curvaton limit, — 1| ~ 2|ry^^ — e*| < 0.1. Note that infiaton mass, ry^^, does not affect the 
spectral tilt in the curvaton limit so the infiaton mass could be of order the Hubble scale without producing a large 
spectral tilt, but it does affect the running of the spectral index. The running (33) in the curvaton limit wx ~ 1 is 
~ as^ which is of the same order of magnitude as the spectral index for 77^0 ~ 1. 
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B. Self- interacting potential 

1. Cosine potential 

We consider an axion-type potential for a weakly broken C/(l)-symmctry (/ ^ Al) [5, 32-34, 40] 



V{x) = 



1 — cos 



(70) 



For X* ^ f effective potential reduces to the quadratic potential (63) with = M / f, but the cosine potential 
has self-interaction terms which become significant for x* ^ f- By symmetry it is enough to consider the range 
^ X*// ^ Ti" for the curvaton VEV during inflation. 

Figure 4 shows the parameter defined in Eq. (29), which determines the contribution of the curvaton to the scalar 
power spectrum for a given value of the inflationary energy scale, iJ*, as a function of x* and m-^/T^ for a cosine 
potential with / = 10^^ GeV. For x* ^ / we recover the previous results for the quadratic potential shown in Figure 1. 
(Note that the y-axis is linear in Figure 4 but logarithmic in Figure 1). For values of > / the higher-order terms in 
the potential reduce the potential gradient and hence slow-down the evolution of x- The curvaton has a larger density 
when it decays than it would have done for the same initial VEV in a quadratic potential. Thus increases and Cc 
decreases relative to the same parameter values in the quadratic potential. In particular this increases the parameter 
range for which the curvaton dominates the primordial power spectrum, ec > e*, relative to the quadratic case. 

In Figure 5 we show the non-Gaussianity parameter, /nli £^nd the tensor-scalar ratio, r^, for different curvaton 
parameter values and a fixed slow-roll parameter, e* = 0.02. Bounds on the tensor-scalar ratio no longer place a lower 
bound on the decay time, m^/T^, as the tensor-scalar ratio becomes small when ec is large for x* ^ ^/ where the 
curvaton VEV is close to the maximum of the cosine potential. Large positive values of /nl also become possible for 
X* ^ T^f due to the non- linear evolution of the curvaton field, even though ~ 1. 

For x* ~ T^f and i?,^. ~ 1 we have from Eq. (46) 



r 5 fg"g 



2 

/2 / ""X 



< , (71) 
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FIG. 5: The plot show contour hnes for the non-Hnear parameter /nl (blue lines) and the tensor-to-scalar ratio (thick red 
lines) as a function of the curvaton parameters, x* £^nd m^/F^, for the cosine curvaton potential, Eq. (70), with / — 10^^ GeV, 
and a fixed value of the inflation slow-roll parameter, e* = 0.02. The black broken line corresponds to tc = 0.02. 



and from Eq. (53) 



3NL 



24 



25 fg"'g^ 



9' 



,g"9 



(72) 



2. Hyperbolic- cosine potential 



We also consider a hyperbolic-cosine potential 



Vix) = 



cosh 



(73) 



For X* ^ f the effective potential reduces to the quadratic potential (63) with = AP / f . Self-interaction terms 
which become significant for x* ^ f ^-nd for x* S> / the curvaton field becomes massive during inflation and evolves 
rapidly to smaller values, hence we will assume X* ^ f our discussion. 

We start by studying the relative contribution of the curvaton to the scalar power spectrum characterised by 
parameter ec defined in Eq. (29). Figure 6 shows Ec as a function of x* and m^/F^ for a hyperbolic-cosine potential 
with / = 10^^ GeV. Again, for x* ^ / "we recover the previous results for the quadratic potential shown in Figure 1. 
For values of x* > / the higher-order terms in the potential increase the potential gradient and speed up the evolution 
of X relative to the quadratic potential. The curvaton has a smaller density when it decays than it would have done and 
thus increases and Cc decreases relative to the same parameter values in the quadratic potential. This decreases the 
parameter range for which the curvaton dominates the primordial power spectrum, relative to the quadratic 

case. 

In Figure 7 we show the non-Gaussianity parameter, /nl, and the tensor-scalar ratio, r^, for different curvaton 
parameter values and a fixed slow-roll parameter, e* = 0.02. As in the case of a quadratic potential bounds on the 
tensor-scalar ratio place a lower bound on the decay time, m^^/F^ if < 16e*. 

Large negative values of /nl are in principle possible due to the non-linear evolution of the curvaton field for X* > /■ 
However, just as in the case of positive /nl for the quadratic potential, extremely large values are not possible for 
finite since Sc becomes small for x* ^ / and hence and /nl ~^ for x*// ~^ -I-cxd. The maximum value 

of /nl (for sufficiently small decay rates, such that ~ 1) occurs when we have e^. e*, i.e, at the boundary of 
curvaton and infiaton limits. 
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FIG. 6: The plot show contours of defined in Eq. 29 as a function of the curvaton parameters, x*i the curvaton VEV, 
and m^/r the dimensionless decay time for the hyperbolic-cosine curvaton potential, Eq. (73), with / = 10^^ GeV. Curvaton 
perturbations dominate the primordial scalar power spectrum for e, S> tc- 

IV. DISCUSSION AND CONCLUSIONS 

In our previous woik [10] we have shown how observables such as the tensor-scalar ratio, r^, and non-linearity 
parameter, /nl, are related to curvaton model parameters, specifically the curvaton VEV, x*, and the dimensionless 
decay rate, F^/m^. In this paper we have allowed for the presence of primordial perturbations due to adiabatic 
inflaton field fluctuations in addition to isocurvature curvaton field fiuctuations during inflation. This introduces 
an additional model parameter, the slow-roll parameter e*, which determines the primordial power spectrum due 
to inflaton field fluctuations relative to the tensor power spectrum. We have constructed an equivalent parameter, 
Ec, which determines the primordial power spectrum due to curvaton field fluctuations relative to the tensor power 
spectrum. For ec <C e* the curvaton fluctuations dominate the primordial scalar power spectrum, ~ 1, and we 
recover the results of our previous work [10]. For 3> e* the inflaton fluctuations dominate the primordial scalar 
power spectrum, ^ 1. 

In practice we have presented two-dimensional contour plots of the tensor-scalar ratio, r^, and non-linearity pa- 
rameter, /nl, as functions of x* and F^/to^ for fixed values of e*. We have shown that a curvaton can produced 
detectable non-Gaussianity and/or gravitational waves for a range of model parameters, even allowing for the presence 
of inflaton perturbations. For a small slow-roll parameter, < Ec, very large values of the non-linearity parameters 
are suppressed (/nl oc (e*/ec)^, ffNL oc (e*/ec)^, etc). Nonetheless /nl may still be produced for < Cc when <C 1 
or in the presence of self-interactions and non-linear curvaton field evolution, \g"g/g''^ \ ^ 1. 

To differentiate between different scenarios for the origin of non-Gaussianity we should examine further the statistics 
of the primordial density field. For example, in the absence of curvaton self-interactions the scale-dependence of the 
non-linearity parameter is given by [13, 37-39] 

Note that if the curvaton dominates both the power spectrum and the higher-order correlators, ~ 1, then /nl is 
independent of scale. If the power spectrum is dominated by inflaton perturbations {w^ ^ 1), such that tiq = n^, then 
the bispectrum and higher-order correlators are still dominated by the curvaton perturbations. Hence we generally 
expect a scale-dependence of the non-linearity parameters /nl and (/nl that since they determine the higher-order 
correlators relative to the power spectrum. The higher-order correlators and the power spectrum inherit different 
scale-dependence from the curvaton and infiaton perturbations respectively. In terms of slow-roll parameters 



'^/nl - 4(1 - w^){2e^ +Vx- V4>) 



(75) 
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FIG. 7: The plot show contour hnes for the non-Hnear parameter /nl (blue lines) and the tensor-to-scalar ratio (thick red 
lines) as a function of the curvaton parameters, x* ''■nd m^/T^, for the hyperbolic-cosine curvaton potential, Eq. (73), with 
/ = 10^® GeV, and a fixed value of the inflation slow-roll parameter, e, = 0.02. The black broken line corresponds to tc = 0.02. 



which is small if we assume slow-roll for both the curvaton and inflaton. However in this case only the inflaton tilt 
is constrained by current observations of the power spectrum and if the curvaton scale-dependence is large then rif^-^ 
could be large. 

The primordial trispectrum also gives important clues about the origin of non-linearity. Figure 8 shows the trispec- 
trum parameters ^nl and tnl as a function of curvaton parameters for a quadratic curvaton potential. In the absence 
of self-interactions non-Gaussianity only becomes large when <^ 1 and in this limit we have 

lOw^X f 36 2 fr.f.\ 

5NL — 7. — JNL , Tnl — 7^3 JNL • ('Dj 

Even allowing for a mixed inflaton+curvaton model with < < 1, we can eliminate to obtain a consistency 
relation between the bispcctrum and trispectrum parameters in this case: 

Snl^nl — — ^/nl^ ■ (77) 
5 

Note that from (76) we can deduce that tnl > O.I^nl^ for a quadratic curvaton potential. If both (/nl and tnl are 
large then the curvaton potential must include self-interaction terms [36]. Such self-interactions can give rise to large 
scale-dependence of /nl and (7nl even in the curvaton-dominated limit, ~ 1 [37-41]. 

Unlike single-field models of inflation, the predictions of the curvaton model are dependent on the initial value of 
the curvaton field during inflation. Although it may not be possible to identify a unique initial value for the curvaton, 
we may be able to specify the expected probability distribution for different models. 

If we take a stochastic approach for the distribution of the curvaton VEV [42] then for a quadratic curvaton 
potential, and assuming that inflation lasted long enough (and assuming a light curvaton, < iJ^), we find a 
Gaussian distribution with variance 

(x^>-^4- (78) 

Huang [43] has argued that in this case, a detection of primordial non-Gaussianity then places a lower bound on the 
tensor-scalar ratio. On the other hand for the self-interacting potential considered here there is an additional scale, 
/, in addition to the effective mass, m^. For a cosine type potential one obtains an almost uniform distribution 
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FIG. 8; Contour lines for the trispectrum non-linear parameters ijnl (thin black lines) and tnl (thick green lines) as a function 
of the curvaton parameters, x* and m^/F^, for the quadratic curvaton potential, Eq. (63), with a fixed value of the inflation 
slow-roll parameter, e* = 0.02. 

for < X* ^ ^/ assuming inflation lasts long enough and is at a high-enough energy scale, while for a hyperbolic 
potential which becomes steep for x* 3> / we expect values with x* ^ / to be suppressed. It would be interesting to 
interpret observations from future observations of primordial non-Gaussianity and/or tensor-scalar ratio in terms of 
curvaton model parameters, incorporating a prior probability distribution for the curvaton VEV to marginalise over 
at least one unknown model parameter. 
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